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The effect of a soft phase core appearance in the center of polytropic star 
is analyzed by means of linear response theory. Approximate formulae for 
the changes of radius, moment of inertia and mass-energy of non-rotating 
configuration with arbitrary adiabatic indices are presented, followed by an 
example evaluation of astrophysical observables. 
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1. Introduction 

In spite of lack of precise Earth-based experimental data involving mat- 
ter at densities higher that the nuclear saturation density the astrophysical 
observations of compact objects give us an unique chance of understanding 
the underlying physics. At present it cannot be excluded that at certain 
density phase transition in nuclear matter will produce a state not observed 
in laboratories e.g. pion or kaon condensation, or de-confined quarks (see 

for review). If the phase transition happens to be of first order, the con- 
sequences from the point of view of observations are extremely interesting. 
Namely, in a first order phase transition, the new phase arise only by nu- 
cleation. This means that the meta-stable core formed during the compact 
star evolution (e.g. accretion, spin-down) nucleates to a stable new-phase 
core. The transition is therefore accompanied by star-quake - the radius 
change, energy release and possible other violent phenomena. 

This work is an extension of the linear response theory developed by 
P. Haensel, L. Zdunik and R. Schaeffer 0, 0] and employed to the case of 
first order phase transition from one pure phase (normal, N-phase) to other 
pure phase (super-dense, S-phase). While one demands that the transition 
should conserve locally electric charge and the number of particles, it occurs 
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at constant pressure, technically by so-called "Maxwell construction" In 
the presence of the gravitational field it is accompanied by a density jump 
at the boundary of phases. 

Due to the works of N. Glendenning 0, it became clear that relaxing 
the condition of local electrical neutrality permits for a coexistence of both 
phases within a range of pressures - a structural mixed-phase transition. 
The volume fraction occupied by the higher-density phase increases from 
zero at the lower pressure boundary to one at the upper pressure boundary. 
In realistic matter, if the surface tension and Coulomb contribution to the 
energy is not to large, the mixed-phase is preferred over a pure phases state. 

Here we derive the formulae suitable for description of the change of 
the parameters (radius, moment of inertia, mass-energy) of a compact star 
under the transition from a meta-stable core to stable core composed of the 
"structural mixed-phase". These changes are proportional to the specific 
powers of the newly-born core. We will neglect the rotation of the star. 
For simplicity, the equation of state (EOS) of matter in both N and the 
mixed phases will be approximated by polytropes. The methods used here 
are similar to those used in a somewhat more complicated case of realistic 
N-phase EOSs presented in 0]. 

The article is arranged as follows: Sect. [21 provides brief description of 
the theoretical background and methods used, Sect. |3] contains results of 
numerical calculations, and Sect. 0] includes conclusion and final remarks. 

2. Linear response theory 

We will describe theoretical background of star linear response to the 
appearance of the soft "mixed-phase" core. The calculation is based on 
expressing the change in the density profile, due to the presence to a small 
core, as the combination of two independent solutions of the linearly per- 
turbed equations of stellar structure [3, 0] • The presence of a denser phase 
in the core changes the boundary condition at the phase transition pressure 
Po (Fig. ^) and allows to determine the numerical coefficients in the expres- 
sion for the density profile change. The leading term in the perturbation of 
the boundary condition at the edge of the new phase results from the mass 
excess due to the lower stiffness, and higher density of the new phase as 
compared to those of the (normal, less-dense) N-phase. 

We assume that at a central pressure P c = P cr it the nucleation of the 
S-phase in a super-compressed core of radius r N , of configuration C, initiates 
the phase transition and formation of the "mixed-phase" core of radius r m 
in a new configuration C*, as shown on Fig. ^ Transition to a mixed 
phase occurring at r m is associated with a substantial drop in the adiabatic 
index of matter, defined as 7 = (n b /P)dP/dn b , from 7 N to 7 m . In realistic 
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Fig. 1. Schematic plot of central pressure P c versus central matter density p c for 
configurations based on a pure N phase EOS and an EOS with a mixed-phase seg- 
ment. The solid line denotes stable states, the dashed line - the states which are 
meta-stable with respect to the transition to a mixed-phase state. For a critical 
central density p cr it the S phase nucleates in the super-compressed core of con- 
figuration C, which results in the transition into a stable configuration C* with a 
mixed-phase core and a central density p*. Configurations C and C* have the same 
baryon number A. 



EOSs, mixed phase is softer than the pure one, because the increase of mean 
density is reached partly via conversion of a less dense N phase into denser 
S phase, and therefore requires less pressure than for a pure phase. This 
remains true for any fraction of the S phase, and leads to a discontinuity 
of 7 at the phases boundary, pq. Realistic examples of such transitions are 
given by UH, ^} and [ill . Table 9.1]. In all those cases, near the boundary 
logarithm of pressure depends linearly on the logarithm of baryon density; 
it clearly indicate that the polytropic approximation of mixed-phase is valid 
in the small mixed-phase core regime. We will henceforth benefit from 
expansion in powers of the (tiny) core radius r m . 

As far as global properties of stars are concerned, the hydrostatic stars' 
equilibria corresponding to EOSs with and without the phase transition 
must be compared. The models are non-rotating, spherically symmetric so- 
lutions of Einstein's equations, usually called Tolman-Oppenheimer-Volkoff 



(TOV) equations [lj, |13| . The moment of inertia was calculated using the 



slow-rotation approximation [lj]. Hydrostatic solutions are labeled by their 
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central density p c . The configurations based on the two EOSs are identical 
up to p c = po; the configuration with such central density will be denoted 
by Co and will be called a "reference configuration". 

As was demonstrated by 0], the leading changes due to a phase tran- 
sition are proportional to the fifth power of r m in the case of stellar radius 
R and the moment of inertia I and seventh power in the case of the grav- 
itational mass-energy E = Mc 2 . During calculations, we assume constant 
7 m in the "mixed-phase" - in principle, the inclusion of r m - dependent 7 m 
contributes to the higher order terms, but this contribution is negligible in 
the case of realistic EOSs (see Q for details). 

When the central density exceeds po, the models begin to differ due 
to the appearance of a softer "mixed-phase" core in configurations corre- 
sponding to the mixed-phase segment in the EOS. For two stars composed 
of equal number of baryons A, greater than a baryon number Aq of the ref- 
erence configuration Co, we compare the global parameters - mass-energy, 
radius, and moment of inertia. Their difference corresponds to the changes 
in these parameters implied by the phase transition in the stellar core. 

The new phase influences the boundary conditions by the presence of the 
prefactor (7 N /7 m — 1), which is qualitatively the same as (p^/Pn — 1) in the 
case of a density jump from pure N phase to pure S phase |2|, La]- Moreover, 
the results obtained by indicate that linear response effects should be 
proportional to the core radius power greater by two than in the case of the 
transition between pure phases. One would expect that the relative change 
of the stellar parameter Q = R, I, E takes the following form: 

SQ = ^ -A3(7 N /7 m " l)(fj, (1) 

where f m = r m /i?o, the power I = 5 corresponds to the radius R and moment 
of inertia I, and I = 7 for the mass-energy E = Mc 2 . The coefficients /3q 
are then the functionals of a reference configuration, i.e. /3q = 0q(Cq). The 
validity of the above statements will be confirmed in the next section by 
means of numerical calculations. 

3. Results 

In order to obtain quantitative results the polytropic EOSs for the N 
and mixed phases will be used. The main reason why we have chosen the 
polytropic EOSs are precision of numerical calculation and simplicity dur- 
ing the exploration of the parameter space. Discussion of the polytropic 
EOSs and their application to relativistic stellar structure calculations was 
presented by Details needed for the calculations are given in the Ap- 
pendix. The Sect, |2] described the form of the expressions for the changes 
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Fig. 2. The response coefficients Pi, Pa, (3e (clashed, solid and dotted line, the 
value of Pe was multiplied by a factor of 10) plotted against the stellar mass for 
a sample polytropic EOS, 7 N = 2.5, 7,„ = 1.5, K N = 0.025 (see the Appendix for 
details). The vertical dotted line denotes the maximal mass of the star, 2.37 M & . 



of stellar parameters obtained in the linear approximation. According to 
Eq. (^P), the coefficients split into two factors, first one depending on the 
mixed phase via 7 S , and the second one (5q being a functional of the ref- 
erence configuration Co- As it is a complicated task to write an analytical 
form of the functional Pq(Cq) we will restrict ourselves to the derivation of 
approximate, but accurate fitting formulae, based on the results of precise 
numerical calculations. 

The Pi, Pr and (3e coefficients obtained in the limit of vanishing mixed- 
phase core r m — > as a function of the mass of the sample reference 
configuration M(Cq) are plotted in Fig. [2J Whole range of masses of the 
reference configuration is presented on account of our ignorance about the 
density at which the phase transition takes place. Results were obtained for 
specific choices of 7 N and r y m . The values of 7 N and K N were chosen in such 
a way, that the EOS of the N-phase produced massive (M > M Q ) neutron- 
star models similar to those obtained for realistic stiff EOS of dense matter. 
In particular, this EOS yields M max = 2.37 M and -RM max = 12.52 km. 
Let us stress, however, that while this EOS is a reasonable representation of 
the EOS of matter with p > 2po, it is completely unrealistic at sub-nuclear 
densities and for masses much smaller than one solar mass. 

One notices a characteristic behavior of (3r and Pi. Despite the fact 
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that we do not know the precise density of onset of the new phase, this 
information is not substantial to predict stars response. The coefficients are 
almost constant for a wide range of masses for /3j and (5r. The behavior of 
Pe is different: for M < 0.8M max this coefficient is proportional to M(Cq). 

The "plateau" values of the coefficients /3j and (5r do not depend on 
the pressure coefficient K N , but on 7 N only. The formulae which describe 
accurately those response parameters read 

0.015 x 7n 9 - 4 0-12 x 7N 9-i 

PRW) - ( 7n 1.13 _ ^8.2 ' PIW) - ( 7n 1. 22 _ !)7.5 ' \ l > 

for /3r and /?/, respectively. 

The case of the mass-energy parameter (3e is different, as it can be seen 
on Fig. |2 As we mentioned before, sufficiently far from the maximal mass, 
it is proportional to M(Cq). We can approximate the value of (3e(Cq) as 
follows: 

a In > a085 ( M \ rtl 

&(°») - v , 6Km °-^-> x Km-J ■ (3) 

It has to be mentioned that the fitting formulae from Eqs. (PHSJ have been 
checked against the adiabatic indices in the range 5/3 to 3.5. The fitted 
expressions are fairly accurate within a wide range of masses, to within 
few per-cent compared to the exact numerical calculations. The response 
coefficients underestimate the magnitude of the linear response for config- 
urations near the maximum allowable mass, M max . The increase of (3q for 
M(Co) — ► M max is due to a "softening" of the reference configuration by 
the effects of General Relativity (similar to the case presented in 0, 0]). 
Thus the linear approximation and, what follows, the approximate expres- 
sions presented above cease to be valid in this region. 

It should be also expected, as far as the realistic EOS of the crust 
is concerned, that the assumption of constant 7 in the crust is an over- 
simplification of the problem. Therefore, the region of small masses (smaller 
than 0.5 Mq) is affected by this unrealistic type of crust. In the case of real- 
istic NSs, the behavior of the [5q coefficients near the star's minimum mass 
should, due to the softening of matter, be generally similar to those near the 
maximal mass (such calculations for the case of realistic EOSs are presented 
in 0). 

The fitting formulae for (3q in the region of validity of the linear-response 
approximation allow us to compute the change of the interesting stellar 
parameters for every pair of the parameters 7 N and 7m . For example, the 
1.4 Mq Co configuration from Fig.[2]has radius Rq = 15.18 km and moment 
of inertia Iq = 2.27 x 10 45 g • cm 2 . The response coefficients are then equal 
(3r = 0.62, (3i = 2.23 and Pe = 0.03. The appearance of r m = 1 km core 
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with 7 m = 1.5 adiabatic index inside the 7 N = 2.5 polytrope star changes 
the radius by ~ 1 cm, a number which can be related to the change of 
radius during macro-glitches in pulsar timing. A relative change of the 
moment of inertia (AI/I = — AO/fi) implies speed-up of the order of 10 -6 , 
again a value common to pulsar macro-glitches. The released energy equals 
3 x 10 44 erg. If the core radius is 4 km (still well described by the linear 
response theory) the change in radius is 8 m, an impressive number even 
in view of the size of terrestrial earthquakes. The value of Afl/fi will be 
equal approximately 10 -3 , three orders of magnitude larger than in biggest 
macro- glitches. Finally, the creation of a 4 km "mixed-phase" core produces 
~ 5 x 10 48 erg of energy. 



The appearance of soft core was studied by means of the linear response 
theory to the star parameters. Simple model of a polytropic star with the 
"mixed-phase" core provides a set of well-approximated formulae, usable for 
estimation of the change of stellar parameters (radius, moment of inertia, 
emitted energy) for a given "stiffness" i.e. the adiabatic index 7 of a star, 
the newly-born core, and its a priori unknown radius. The fitted formulae 
are precise to within a few per-cent. Resulting parameter changes are of 
the order of observed astrophysical phenomena, which gives hope for future 
observations of phase transitions inside real neutron stars. Inclusion of a 
varying adiabatic index certainly modifies the values of the response coef- 
ficients, presented in 0, but the comparison with the results shown here 
proves that the difference is not dramatic. 
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The relativistic polytrope is defined as the power-law dependence be- 
tween pressure P and baryon number density n b (see, e.g. |l5j): 



where 7 is often called adiabatic index, and K is the pressure coefficient . 

1 The coefficient K is expressed in pc 2 /n 7 units, where p := 1.66 x 10 14 g/cm 3 , and 
n := 0.1 fm -3 . 



4. Conclusions 
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Assume that matter is strongly degenerate, so that T = approximation 
is valid. For simplicity, we consider one type of baryons in the outer N phase, 
with the rest mass m N = 1.66 x 10 -24 g. The total mass-energy density £ 
of particles of rest mass m N in N phase is related to their baryon number 
density n b by the First Law of Thermodynamics 

K 

£(n h ) = — nl N + m N c 2 n b . (5) 

7n — 1 

Baryon chemical potential match to the change of the energy of matter, 
at constant P and T = 0, due to an increase of the baryon number by one. 
This implies 

d£ P + £ K N 7 N „ i 2 

M n b) = — = = 7< + "InC • (6) 

dn h n b 7 N - 1 

At zero pressure (the surface of the star) the chemical potential \i is equiv- 
alent to particle rest energy. 

The connection between two polytropes (the one corresponding to the 
core being softer than the outer part, 7 m < 7 N , see Fig. ^and discussion) 
aims for approximating the transition to the mixed-phase, known from many 
realistic calculations 0, Hoi . UlT ] . Thermodynamic equilibrium implies that 
the choice of parameters K m , 7 m , and the mean mass m m should assure the 
continuity of pressure and baryon chemical potential are continuous along 
the transition point n , that is 

TV" _ T/ 7N~7m _ _ WN — 7m)P(^o) ,„s 

(7n - l)(7m - l)n c 2 
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